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We theoretically investigate the doping evolution of the electronic state of high-Tc cuprate 
on both sides of the half-filling on the basis of the three-dimensional three-band Hubbard 
model with a layered structure using the Hartree-Fock approximation. Once a small amount 
of holes or electrons are doped into the half-filled state, our model exhibits the charge-transfer 
insulator-to-metal transition along with a chemical potential jump. At the same time, the 
doped holes or electrons are inhomogeneously distributed, and they tend to form clusters in 
the vicinity of the half-filling. This suggests the possibility of microscopic phase separation 
with the separation between the metallic and the insulating regions. 


1. Introduction 

As shown by the photoemission spectroscopy, the undoped high-Tc superconducting 
cuprate (HTSC) is a charge-transfer insulator, where the Cu 3d electrons are almost lo¬ 
calized by the strong electron correlation.^ When the rare-earth element is substituted out 
of the two-dimensional (2D) Cu02 layers and the number of holes or electrons doped into 
the Cu02 layers is increased, Cu 3d electrons hybridized with O 2p electrons achieve itin¬ 
erancy and display superconductivity. Moreover, in slightly hole-doped La 2 -a;Srj;Cu 04 with 
0 < X < 0.12, the chemical potential shift suppression is observed by photoemission spec¬ 
troscopy (PES).^’^ The electronic phase separation between the antiferromagnetic insulating 
phase and the superconducting phase is considered to be one of the reason for the suppres¬ 
sion of the shifting of the chemical potential.'^ The phase separation assumes that the doped 
carriers are inhomogeneously distributed due to the strong electron correlation. 

Many experimental findings have suggested that the electrons under such circumstances 
favor some types of spontaneous ordering in certain doped regions. For instance, in or¬ 
der to elucidate the anomalous suppression of the superconducting transition tempera¬ 
ture in Lai. 875 Bao.i 25 Cu 04 ,^ the spin and charge correlations in Lai, 875 Bao.i 25 Cu 04 or 
(La,Nd) 2 - 3 ;Srj;Cu 04 with x « 0.125 have been intensively studied, and the stripe order has 
been observed by neutron scattering.The stripe order has also been observed by x-ray 
scattering.Furthermore, an electron paramagnetic resonance study showed that micro- 
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scopic electronic phase separation occurs in La 2 -xSr 2 ;Cuo. 98 Mno.o 204 with 0.01 < x < 0.06/® 
and a Cu nuclear magnetic resonance (NMR) study suggested that a large charge droplet 
(’blob’) is formed in the electron-doped Ndi. 85 Ceo.i 5 Cu 04 _ 5 .^^ 

Much theoretical works has also been performed to study the behavior of the doped carri¬ 
ers in HTSC. Pioneering works adopting the Hartree-Fock approximation (HFA) have studied 
the stripe order in Lai.875Bao.i25Cu04 on the basis of the 2D one-band Hubbard model^^’^^ 
or the 2D two-band Hubbard model.Furthermore, dynamical mean field theory (DMFT) 
has been exploited to study the stripe phase on the basis of the 2D Hubbard model with L 
non-equivalent sites, where L = 8,..., 160.^® The DMFT approach has also been adopted to 
analyze the three-band Hubbard model.In some of these works,the DMFT approach 
was combined with the local density approximation (LDA). Another study^® considered the 
possibility of two-sublattice antiferromagnetism. All of these works have successfully repro¬ 
duced the Zhang-Rice singlet band.^^^^^ This shows that the three-band Hubbard model is an 
appropriate model for HTSC near the half-filling and that the DMFT is a powerful tool for 
analyzing its electronic state. However, when we investigate the inhomogeneous electronic dis¬ 
tribution near the half-filling, we need to adopt the model for a large number of non-equivalent 
sites. In general, the DMFT costs much more than the HFA to analyze the model with a large 
number of non-equivalent sites. 

In this paper, we analyze the normal ground state of the 3D three-band Hubbard model 
with a single-layered perovskite structure in order to study the evolution of the electronic state 
when holes or electrons are doped into the undoped HTSC. We consider 256 non-equivalent 
copper sites for each rectangular parallelepiped super cell, and adopt the HFA for these con¬ 
ditions. We performed the calculation without any assumptions about the electronic distribu¬ 
tion, and we obtained fully self-consistent solutions except near the 1/8-filling. These solutions 
showed the chemical potential jump at half-filling, which means that the electron suddenly 
becomes itinerant when a small number of holes or electrons are doped. Moreover, the doped 
holes or electrons tend to form clusters in the vicinity of the half-filling. These clusters are con¬ 
sidered to form a metallic region, and are surrounded by the insulating region. This suggests 
the possibility of microscopic electronic phase separation in HTSC near half-fillig. 

2. Formulation 

Our 3D three-band Hubbard model Hamiltonian, H, is composed of d-electrons at each 
Cu site and p-electrons at each O site. To consider the spatial inhomogeneity, we introduce 
the rectangular parallelepiped super cell containing Nc Cu and 2Nc O sites as a unit cell. 
Thus, H is defined as follows; 

Afc iVc .. Afc iVc 

A^ko-'^b'kA'kcr + vEE '^lk+qt'^Ik'-q4,'^*kq'ijkt EEcTOk,. (1) 

i=l j=l k(T i=l kk'q i=l kcr 
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Here we use the abbreviations = {dlkcr Pika P^L) = Xdika Pik^ Pika)^ where 

dikaid}]^^) and are the annihilation (creation) operators for the d-orbital and 

p*(j/)-orbital electron on the i-th site, as specified by the momentum k and spin a = {tj-i 
}, respectively. U, N, and fj, are the on-site Conlomb repnlsion between the d-orbitals, the 
number of k-space lattice points in the first Brillouin zone (FEZ), and the chemical potential, 
respectively. The FEZ is defined in the reciprocal space to the lattice whose unit cell contains 
Nc Cu and 2Nc O sites. The unit cell is schematically shown in Fig. 1. The two non-equivalent 
Cu02 layers, indicated with LI and L2, are alternatively stacked along z-axis. On each Cu02 
layer, the size of the unit cell along x-axis and y-axis is and Ny Cu sites, respectively. 
Thus, when we set = 8 and Ny = 16, = 2 x 8 x 16 = 256. We take the primitive 

translation vectors for the unit cell to be ±8a(x-|-y), ±8a(x —y), and ±2cz, where x, y, and 
z are the unit vectors for the x-axis, y-axis, and z-axis, respectively. in Eq- (1 ) is defined 
as follows: 
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where A^p is the hybridization gap energy between the d- and -orbitals. where p G 
(d, p*, p^), is the coordinate of the p-orbital electron on the f-th site. is Kronecker’s 

delta, i.e., it is 1 for = R^y,' and 0 for R*^ ^ Rjy>'. In the following, we take both a and 
c to be a unit of length and set a = c = 1. Then, we can represent 
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where tpp is the transfer energy between a p*-orbital and a p^-orbital, t^p is that between a 
d-orbital and a -oxhiisl, and t±_ is that between d-orbitals, respectively. In this study, tdp 
is the unit of energy. 
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Fig. 1. The schematic figure of the unit cell. LI and L2 indicate the two non-equivalent Cu02 layers 
alternatively stacked along z-axis. On each Cu02 layer, the size of the unit cell along x-axis and 
y-axis is and Ny Cu sites, respectively. In total, the unit cell contains Nc Cu and 2Nc O sites. 


where Nc = 2 x x Ny. 


We adopt the HFA with respect to every Nc Cu site and two spin states, and we only 
consider collinear spin states. Thus, we define 


and we approximate Eq. (1) as follows: 
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Then, we conduct a self-consistent calculation of Eqs. (3), (4), and (5) and obtain self- 
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consistent fields riida, and Uipya, where we define 
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3. Results and Discussion 

In the numerical calculations, we divide the FEZ into N = 16 x 32 x 4 equally-spaced 
rectangular parallelepiped. The parameter sets are selected as tdp = 1.0eV, tpp = —0.3eV, 
t± = 0.005 eV, Adp = 0.0 eV, and U = 6.0 eV. In order to obtain self-consistent solutions, we 
need to set initial values for fields Uida, n-ip^a, and Uipya and carry out iterative calculation until 
all of these fields have sufficient accuracy. Here, we chose the initial values for the fields from 
the uniform-distributed random numbers. In this manner, we obtained the solutions for the 
doping region near half-filling; 0.0 < (5h E 0.1 and 0.0 < (5e < 0.1, where 5e = —(Jh- For these 
solutions, all of the obtained fields riida, n-ip^a, and Uipya- have three digits of accuracy. Once 
they get three digits of accuracy, all of these fields rapidly converge. Thus, we can consider 
the solutions with these fields as fully self-consistent. We also tried to obtain the solutions 
for other doping region e.g. (5h ~ 0.125 but failed. For these doping regions, the ground state 
accompanied by certain types of long-period superlattice structure would be stable. The fields 
corresponding to such ground state could be hardly obtained by our iterative calculation. We 
henceforth concentrate our discussion on the doping region near half-filling. 

The doping dependence of the chemical potential ^ for our fully self-consistent solutions is 
shown in Fig. 2. The chemical potential changes rapidly at half-filling, which is consistent with 
experimental results obtained by comprehensive PES studies on HTSC.^® The magnitude of 
this change is about 3.0 eV, and the value is in accordance with Uf2 — Adp, which is equal to 
[7/2 in our calculations. This fact can be explained by the doping dependence of the density 
of states (DOS), which is shown in Fig. 3. In this figure, we show the mean DOS over the 
solutions with similar doping and chemical potential as DOS for each doping state, and we 
indicate each doping state with the mean doping over these solutions. All these labels in Fig. 3 
are summarized in Table I For instance, DOS for 5e = 0.094 in Fig. 3 is the mean DOS over the 
seven solutions for which we found 5e = 0.094 ± 0.000 and ^ = 6.301 ± 0.015 as a mean with 
standard error of the mean of doping and chemical potential, respectively. In the slightly 
electron-doped case, Se = 0.094,0.067,0.031, the Fermi level crosses the upper Hubbard ri¬ 
band, while in the slightly hole-doped case, rih = 0.032,0.066,0.095, the Fermi level crosses 
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Fig. 2. (Color online) The doping dependence of the chemical potential /i. (5e and dh are the numbers 
of doped electrons and holes per Cu 02 unit, respectively. 




Energy (eV) 

Fig. 3. (Color) The doping dependence of the mean DOS. The Fermi level at zero energy is indicated 
by the vertical dashed line. The lines of upper half and lower half indicate the DOS for spin up 
and spin down, respectively. In the ascending orders of energy, we attribute each three blocks of 
DOS to p-band plus lower Hubbard d-band, p-band, and upper Hubbard d-band, respectively. 

the p-band. The upper d-level is about Uriidal'^ larger than the p-level when A^p = 0.0 eV, 
due to the Hubbard splitting. In the energy range between the p-band and the upper Hubbard 
d-band, the density of states is almost zero. In the undoped case, dh = 0.000, the Fermi level 
locates in this energy gap, and it causes the chemical potential to jump by U/2 at half-filling. 
The quantum Monte-Carlo calculation of the two-band Hubbard model in the limit of infinite 
dimensions gave the same result for the chemical potential jump as our calculation.^® Thus, 
the chemical potential jump at half-filling is the characteristic behavior of the multi-band 
Hubbard model composed of both Cu 3d electrons and O 2p electrons independent of the 
dimensionality of the model. 
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N 

5(e,h) ± SEM 

^±SEM 


= 0.094 

7 

0.094 ± 0.000 

6.301 ±0.015 

<5e 

= 0.067 

8 

0.067 ±0.000 

6.287 ±0.011 

-5e 

= 0.031 

8 

0.031 ± 0.000 

6.246 ± 0.006 

4 

= 0.000 

5 

0.000 ± 0.000 

4.823 ±0.019 

4 

= 0.032 

8 

0.032 ± 0.000 

3.054 ± 0.005 

4 

= 0.066 

8 

0.066 ± 0.000 

2.971 ±0.011 

4 

= 0.095 

8 

0.095 ± 0.000 

2.892 ± 0.004 


Table 1. The statistical values of the solutions with similar doping (5(e,h) chemical potential fx, 
which are selected from the ones in Fig. 2. N and SEM mean number of solutions and standard 
error of the mean, respectively. These labels are commonly used in Figs. 3, 4, 5, 6, and 7. 


The inhomogeneous distribution of every obtained field rii^a can be observed in our so¬ 
lutions. The doping dependence of the distribution of the obtained fields can be shown by 
histogram, as in Fig. 4. In this figure, we show the accumulated histograms over the solutions 
with similar doping and chemical potential as the histogram for each doping state, and we 
indicate each doping state with the mean doping over these solutions as well as in Fig. 3. At 
half-filling, where the ground state is insulating, almost all riida have the same value near 1 
and almost all riipx^ and njps/o- have the same value near 2. In the electron-doped case, with the 
increase of the doped electrons, only the peak for riida broadens and its center shifts higher. 
In contrast, in the hole-doped case, with increasing hole density, not only the peak for riida 
but also the peaks for riipx^ and njpHo- broaden and their centers shift lower. Except for their 
variances, the doping dependence of the average of riida is similar to that of Ud obtained by 
the LDA-I-DMFT calculation, and the doping dependence of the average of or njpHo- is 
similar to that of npP'^ 

The spacial distribution of the magnetic moment rriip = also becomes inhomo¬ 

geneous by doping. In the same manner as in Fig. 4, the doping dependence of the distribution 
of the magnetic momenta magnitude can be shown by histogram in Fig. 5. At half-filling, one 
half of niid is exactly at 1 and the other half is exactly at —1. This indicates that the d-electron 
spins are fully polarized. Both in the electron-doped case and in the hole-doped case, with an 
increase of the doped carrier, the two peaks for niid become broader and their centers shift 
toward zero. This means that the number of the fully-polarized d-electron spins decreases and 
that the number of the doubly-occupied Cu sites increases. We should note that the tails of 
rriid in the electron-doped case extend more than those in the hole-doped case. 

The difference between the electron-doped case and the hole-doped case is caused by the 
doped carriers being differently distributed in the unit cell. In order to show this distribution, 
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Fig. 4. (Color) The doping dependence of the histogram for the number of electrons per site. We can 
attribute the peaks whose number of electrons per site are at or near 1 to the ones for riida and 
those whose number of electrons per site are at or near 2 to the ones for riipx^ or rupya- 



Fig. 5. (Color) The doping dependence of the histogram for the magnetic moment per site. We can 
attribute the peaks whose magnetic moment per site are at or near ±1 to the ones for rriid and 
those whose magnetic moment per site are at or near 0 to the ones for rriip^ or rriipv. 


we define the distribution function of the doped carriers as follows: 


4(r) = 


(1 




(9) 


The doping dependences of 5h(r) for the electron-doped case and for the hole-doped cases are 
shown in Figs. 6 and 7, respectively. In these figures, the red spots indicate the location of 
doped holes and the blue ones indicate the location of doped electrons. In the electron-doped 
case, as shown in Fig. 6, the doped electrons form blobs even in the slightly electron-doped 
case, with 5e = 0.031. These blobs can be identihed as the ones in Ndi.85Ceo,i5Cu04_5, 


8/14 


















J. Phys. Soc. Jpn. Full Paper 

suggested by the Cu NMR study.^® On the other hand, as shown in Fig. 7, the doped holes 
stay within a single CUO 4 cluster in the slightly hole-doped case, with (5h = 0.032. Hence it 
can be recognized that the cluster formed by the doped electrons is larger than the one formed 
by the doped holes. We can quantitatively clarify this tendency with the doping dependence 
of the number of the sites having more than certain value of excess carriers in Fig. 8 since 
such a number should increase with the size of the cluster formed by the doped carriers. 
The difference between the electron-doped case and the hole-doped case is explained in the 
following; In both cases, the doped carriers basically tend to form an extended cluster, since 
the strong on-site Coulomb repulsion hinders the double occupancy on each site, and instead 
align them next to each other to gain more kinetic energy. The study on Coulomb gas ordering 
in a 3D layered system by the Brownian dynamics approach supports our explanation.^^ In 
our lattice model, the size of such a cluster depends on the number of adjacent orbitals where 
the doped carriers are allowed to occupy. In the electron-doped case, the doped electrons have 
room to sit only on the Cu sites, because all O sites are fully filled by electrons. On the other 
hand, in the hole-doped case, the doped holes have room to sit on both Cu sites and O sites. 
Therefore, when the number of doped electrons is as many as the number of doped holes, the 
cluster formed by the doped electrons is larger than the one formed by the doped holes. 

4. Conclusion 

In this paper, we conducted HFA calculations for the 3D three-band Hubbard model with 
a single-layered perovskite structure considering a large number of non-equivalent sites. We 
obtained the fully self-consistent solutions both for the electron-doped and hole-doped cases at 
or near half-filling. Our solutions show the chemical potential jump at half-filling. The jump 
can be explained by the DOS dependence on doping, which is characteristic to the multi¬ 
band Hubbard model composed of both Cu 3d electrons and O 2 p electrons independent of 
the dimensionality of the model. Inhomogeneous electronic distributions near half-filling are 
observed in our solutions. There is a remarkable difference in the inhomogeneous electronic 
distributions between the electron-doped and hole-doped cases. That is, the clusters formed 
by doped carriers extend more in the electron-doped cases than in the hole-doped cases. The 
difference between the electron-doped and hole-doped case is caused by the difference in the 
species of orbitals the electron and hole are allowed to occupy, which should be explained only 
on the basis of the multi-band Hubbard model. Thus, the theoretical approach on the basis 
of the multi-band Hubbard model can explain both the chemical potential jump at half-filling 
and inhomogeneous electronic distributions near half-filling in a comprehensive way. 
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Fig. 6. (Color online) The doping dependence of (5h(r) for the electron-doped case. Each figure is a 
single snapshot of the solutions attributed to each doping state, indicated with the mean doping 
over these solutions. 10/14 
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Fig. 7. (Color online) The doping dependence of (5h(r) for the hole-doped case. Each figure is a single 
snapshot of the solutions attributed to each doping state, indicated with the mean doping over 
these solutions. 11/14 
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Fig. 8. (Color online) The doping dependence of the number of the sites having more than 0.1 excess 
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